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Abstract 

We derive upper and lower bounds on the critical temperature T c and the energy gap S (at zero 
temperature) for the BCS gap equation, describing spin 1/2 fermions interacting via a local two- 
body interaction potential XV (x). At weak coupling A <C 1 and under appropriate assumptions on 
V(x), our bounds show that T c ~ Aexp(— B/X) and S ~ C exp(—B/X) for some explicit coefficients 
A, B and C depending on the interaction V(x) and the chemical potential /i. The ratio A/C turns 
out to be a universal constant, independent of both V{x) and /i. Our analysis is valid for any ^; 
for small /x, or low density, our formulas reduce to well-known expressions involving the scattering 
length of V(x). 
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I. INTRODUCTION 



The recent advances in trapping and cooling of cold atoms have led to renewed interest 
in the behavior of ultra-cold fermionic gases. Under the assumption that the interactions 
among the individual fermions are weak, the system shows a superfluid behavior at low tem- 
perature that is well described [l, 2, 3], 4], 5), [^,3, Is] by the Bardeen-Cooper-Schrieffer (BCS) 
model 9j. The BCS model was originally introduced as a model for electrons displaying 
superconductivity, and has played a prominent role in condensed matter physics in the fifty 
years since its introduction. We shall not be concerned here with a mathematical justifica- 
tion of the approximations leading to the BCS model, but rather with an investigation of 
its precise predictions. 

In this paper, we study the BCS gap equation for fermionic systems with general local 
interaction potentials. We give a rigorous derivation of expressions for the critical temper- 
ature T c and the energy gap H (at zero temperature) that are valid to second order Born 
approximation. More precisely, for all interaction potentials V(x) that create a negative 
energy bound state of the effective potential on the Fermi sphere (see Eq. f[TTj) below; a 
sufficient condition for this property is that L 3 V(x)dx < 0), we show that 

Tc = ^l e -/i^) (1 ) 



where /1 > is the chemical potential, 7 m 0.577 denotes Euler's constant, and < is 
an effective scattering length. Units are chosen such that h = ks = 2m = 1, where m is 
the mass of the fermions. To first order in the Born approximation, is related to the 
scattering amplitude of particles with momenta on the Fermi sphere, but to second order 
the expression is more complicated. The precise formula is given in Eq. ffl6|) below. For 
interaction potentials that decay fast enough at large distances, we shall show that b^ reduces 
to the scattering length a of the interaction potential in the low density limit, i.e., for small 
/1. 

We emphasize that not only are our results mathematically rigorous, but our analysis 
holds for arbitrary (positive) values of the chemical potential fi. In particular, the formula 
(P) generalizes previously known results valid only for small fi, i.e., low density. To the 
best of our knowledge, the correct expression for b^ does not seem to have appeared in the 
literature before. We use precise spectral analysis of the (linearized) BCS gap equation in 
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order to derive our results. 

For simplicity, we restrict our analysis to local interaction potentials, which is the case 
of interest when describing dilute Fermi gases. Our methods are applicable in a much more 
general setting, however, and generalizations to non-local potentials (as used in the theory 
of superconductivity) are straightforward. 

For interaction potentials that have non-positive Fourier transform (ensuring, in par- 
ticular, that the BCS pair wavefunction is unique and has zero angular momentum), 
we shall prove similar results for the zero temperature energy gap, which we denote by 



H = min p a/ (p 2 — /i) 2 + |A(p)| 2 . It turns out that, at least up to second order Born approx- 
imation, 

in this case. This equality is valid for any density, i.e., for any value of the chemical potential 
jj,. In particular, H has exactly the same exponential dependence on the interaction potential, 
described by b^, as the critical temperature T c . 

Despite the huge physics literature concerning the BCS gap equation, rigorous results 
concerning its prediction are sparse. Our analysis here relies on the previous studies in [h]] 



in this direction, and extends the results of 
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Before giving the mathematically precise statements in Section IHIl we shall give a non- 
technical discussion of our main results in the next section. Sections IIVHVII contain their 
proofs. 



II. DISCUSSION OF MAIN RESULTS 



We consider a gas of spin 1/2 fermions at temperature T > and chemical potential 
H > 0, interacting via a local two-body interaction potential of the form 2XV(x). Here, 
A > is a coupling parameter, and the factor 2 is introduced for convenience. The superffuid 
phase of the system is described by the BCS gap equation 

where E{q) = a/ {p 2 — /i) 2 + |A(g)| 2 . The Fourier transform of V(x) is denoted by V(p). 

We are interested in the critical temperature T c which, in mathematical terms, is defined 
by the property that ([3]) has a non-trivial (that is, not identically vanishing) solution for 
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T < T c , while there is no solution for T >T C . In the limit A — > 0, we shall show in Theorem [T] 
that 

lim ( In ( £ ) + ^— T l = 2 - 7 - ln(8/7r) . (4) 

Here, 7 ~ 0.577 denotes Euler's constant. For small A, T c is thus given by (JTJ). The quantity 
6 M (A) plays the role of an effective scattering length. In fact, as \i — > 0, it reduces exactly to 
the scattering length of 2W{x). For general /x, the expression is slightly more complicated, 
however. The precise formula is given in the next section, see Eq. fll6p . In certain special 
cases, including all potentials with non-positive Fourier transform, 6 M (A) is given by 



where </?(x) = (2 / K 2 fi)~' L ' 2 V(x) sm(^yji\x\)/\x\ and </?(p) is its Fourier transform. 

For A the solution of the gap equation (J3j) at zero temperature T = 0, the energy gap 
of the system is given by S = min p a/ (p 2 — /i) 2 + |AQo)| 2 . For small A, it shows the same 
behavior as the critical temperature. More precisely, we shall show in Theorem [2] that 

lim fin (£) + -— ^— — ) = 2 - ln(8) . (6) 



a^o V \EJ 2^/Ib^X) 
Together with this shows the validity of (J2J) for small A. 

We emphasize that our results hold for a large class of interaction potentials V(x) and are 
hence suitable to describe a wide range of physical situations. It is not necessary to make 
the approximation of zero-range, for instance, as is often done in the physics literature. Our 
results can therefore be interpreted as an a posteriori justification of such an approxima- 
tion. Note that V(x) does not necessarily have to be interpreted as the exact interaction 
potential among the particles but can include effective interactions arising from higher order 
contributions not taken into account in the BCS approximation. 

Formulas (CQ) and (j2j) are well-known for small /1, with 6 M (A) replaced by the scattering 
length of 2XV(x), which we shall denote by Oq(A). Our analysis shows that they are valid 
for all n > 0. For given potential V(x), the effect of non-zero /i on 6 M (A) can easily be 
calculated. For instance, for a Gaussian potential V(x) = — Kexp(— k|x| 2 ), the difference 
between 6 M (A) and the scattering length a (A) of 2XV is given by 



MA)-ao(A) 



.3/2 



(^X + V^ ln(l + V2) + ^= A 2 + 0(A 3 )^) for//<K 



Higher order corrections can be calculated from (J5J) as well. 



III. PRECISE STATEMENT OF RESULTS 



In the following, we shall assume that V is a real-valued potential that has some mild 
regularity properties, namely V G L 1 (IR 3 ) n L 3 / 2 (R 3 ), i.e., J R3 \V(x)\ p dx < oo for 1 < 
p < 3/2. In the BCS approximation, the system under consideration is described by the 



BCS functional JF, which was introduced by Leggett in [1]. If 7 denotes the momentum 



distribution of the fermions, and a is the pair wavefunction, JF is given by [10] 

JF(7, a) = / (p 2 — fi) / y(p)dp + A / \a(x)\ 2 V(x)dx — TS(j, a) . (7) 

Here, S denotes the entropy functional 

5(7, a) = - / [s(p) \ns(p) + (l - s(p)) In (l - s(p))] dp, 
Jm. 3 

with s(p) determined by s(l — s) = 7(1 — 7) — |d| 2 . Here and in the following, we use a caret 
to denote Fourier transform; i.e., f(p) = (27r)~ 3//2 J r3 f(x)e~ tpx dx. The BCS functional T is 
naturally defined for pairs of functions (7, a) with 7 G L 1 (IR 3 , (1 +p 2 )dp), < 7(p) < 1, and 
a G -f/^IR^Gfe) satisfying |a(p)| 2 < 7(p)(l — lip))- 

The existence of a minimizer of T was shown in [h]]. The corresponding variational 
equation satisfied by a minimizer can be formulated as follows. Given a minimizing pair 
(7, a), one defines A as 



Hp) = t^tt«(p) • (8) 

2 - 7(P) 



It satisfies the iJCS 1 ^ap equation 



with -E(p) = a/ (p 2 — aO 2 + |A(p)| 2 . If A does not vanish identically (or, equivalently, the 
minimizing a does not vanish identically), the system is said to be in a superfluid phase. 



A. Critical Temperature 



It was shown in [lCj that there exists a critical temperature T c below which the gap 
equation ([9]) has a non-trivial solution (i.e., A does not vanish identically), and above which 



it does not. Equivalent ly, the a minimizing the BCS functional T is identically zero for 



id, 



T > T c , while it is non-zero for T < T c . This critical temperature is characterized by 
Thms. 1 and 2] 

T c = inf {T > : inf spec (K T)ll + XV) > 0} . (10) 
Here, i£i> is the multiplication operator in momentum space 

e (p 2 -H)/T + I 



Kt,, = {p 2 - /*) 



e (p 2 -»)/T _ I 



Note that K?^ > 2T > 0, and lim^o^r,^ — \p 2 — Observe that Eq. (fTUj) characterizes 
the critical temperature for the nonlinear BCS equation in terms of the spectrum of linear 
operators. 

We assume that fi > henceforth. For small coupling A, the critical temperature is 
determined by the behavior of the potential V(x) on the Fermi sphere the sphere in 
momentum space with radius JJL. We denote the uniform (Lebesgue) measure on by 
dui. 

Let V p : L 2 (f2 M ) — > L 2 (f2 M ) be the self-adjoint operator 



0» (P) = (2^^ / V(P " g)M(g) ' (U) 



Since V G L^IR 3 ) by assumption, V"(p) is a bounded continuous function, and hence V M is a 
Hilbert-Schmidt operator. It is, in fact, trace class, and its trace equals ^ L 3 V(x)dx. Let 

e M = inf spec V M 

denote the infimum of the spectrum of V M . Since is compact, we have < 0. We 
note that operators of the form fill!) have appeared in related works on bound states and 



scatte ring properties of pseudo-differential operators, see [12l. Il3|. 



In 



111 ], it was shown that when e M < the asymptotic behavior of T c as A tends to zero 



is, to leading order, given by 



limAln^ =-— . (12) 



That is, T c ~ exp(— l/A|e M |) for small A. In the following, we shall derive the second order 
correction, i.e., we will compute the constant in front of the exponentially small term in T c . 
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Hw» = / d\ P \ 



For this purpose, we define an operator W M on L 2 {VL^) via its quadratic form 

ir^-iiy ^(i^)l 2 -l^p/bi)l 2 ) 

\\p\ 2 - H Us 2 

+ j dSl\^p/\p\)\^j . (13) 



Here, <p(jp) = (27r)~ 3 / 2 J q V(p — q)u(q)duj(q), and (|p|, Q) G K+ x § 2 denote spherical coordi- 
nates for p G M 3 . We note that since V G L X (IR 3 ), J g2 dfi |y3(p)| 2 is Lipschitz continuous in \p\ 
for any u G L 2 (IR 3 ). (See Eq. fl34l) below.) Hence the radial integral in ffl3l) is well-defined, 
even in the vicinity of \p\ ~ JJL. For large the integral converges because V G L 3 / 2 (R 3 ). 
We shall, in fact, see that the operator is of Hilbert-Schmidt class (see the proof of 
Theorem [1] in Section ITVT) . 
For A > 0, let 

and let 6 M (A) denote its ground state energy, 

6 M (A) = inf specif. (15) 

We note that if e M < 0, then also 6 M (A) < for small A. In fact, if the eigenfunction 
corresponding to the lowest eigenvalue e M of V M is unique and given by u G L 2 (f2^), then 

&„(A) = <r*|B» + 0(A 3 ) = A^| - A 2 ^>1 + 0(A 3 ) . (16) 

In the degenerate case, this formula holds if one chooses u to be the eigenfunction of V M that 
yields the largest value (m|W m |m) among all such (normalized) eigenfunctions. 

With the aid of 6 M (A), we can now state our first main result concerning the asymptotic 
behavior of the critical temperature T c for small A. 

THEOREM 1 (Critical Temperature). Let V G L^M 3 ) n L 3/2 (M 3 ) and let /i > 0. 

Assume that e^ = infspecV^ < ; and let 6^(A) be defined in l[T5\) . Then the critical 
temperature T c for the BCS equation, given in Eq. < f77J)) . is strictly positive and satisfies 

lim (in (£) + - J u ^ ] = 2 - 7 - ln(8/7r) . (17) 



TJ 2 y /fib ll {X) 



Here, 7 ~ 0.577 denotes Euler's constant. 
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The Theorem says that, for small A, 

T c ~ ^/(V^mM) . (18) 
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Eq. (fT6|) can be reformulated as 



7T 1 (ttlW^Itt) 



+ x ' 7 +Q(A) (19) 



2^u6 M (A) Ae M y/Ie^ 

as A — > 0, where w is an eigenf unction of V M in ffTTl) with eigenvalue e M < 0. We note that for 
radial potentials V(x), the eigenfunction u corresponding to the lowest eigenvalue e M of V M 
will be an eigenfunction of the angular momentum. It need not have zero angular momentum, 
however, but can, in principle, have arbitrarily high angular momentum depending on the 
details of V(x). 

Because of ffT9l . Theorem [1] could alternatively be formulated as 



It is natural, however, to absorb the last term into the definition of b^(X), as we do here, 
since b^(X) can be interpreted as a (renormalized) effective scattering length of 2XV(x) (in 
second order Born approximation) for particles with momenta on the Fermi sphere. In 
fact, if V is radial and J R3 V(x)dx < 0, it is not difficult to see that for small enough /i 

owest eigenvalue of V M is the constant 



the (unique) eigenfunction corresponding to the 

function u(p) = (47r/i) -1 / 2 . (See Section 2.1 in Uj.) For this u, we have 



\im(u\Bu\u) = (A/4tt) / 

Jk- 



V{x)dx - (A/47r) a / V ^ V ^ dxdy = a (A) . 

\x — y\ 



Here, ao(A) equals the scattering length of 2AV^ in second order Born approximation. As- 
suming a certain decay rate of V at infinity we can, in fact, estimate the difference between 
(u\Bfj,\u) and ao(A). 

PROPOSITION 1 (Relation to Scattering Length). Let V G L^M 3 ) n L 3/2 (M 3 ) ; 
ji > 0, and let u(p) = (47r/i) -1 / 2 be the constant function on the sphere Q^. Let < (3 < 2, 
and assume that J R3 |y(x)||x|^da; < oo. Then 

lim — — ^((ulB^u) - a (A)) = . 



A'- 
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In particular, if (3 > 1, this implies that 

lim — f , , I , ; L_]=o. (20) 

^o^\(u\B^\u) ao(X)J 

As explained above, for radial potentials V with f R3 V(x)dx < the eigenfunction corre- 
sponding to the lowest eigenvalue of V M in ffTTj) is the constant function u(p) = {A-Kjj)^ 1 / 2 
for small /i, and thus 6 M (A) = (ulB^u) + 0(A 3 ) for small enough n in this case. Hence ( IT7|) 
holds, with 6^(A) replaced by (u|£>^|u). If, in addition, f R3 |V(x)||x|(ia; < oo, one can use 
f l20|) to replace by ao(A), and arrive at the approximation 

T c ~ ^ e 7r/(2v^a (A)) 

C 7T 

for the critical temperature for small A and small fi. This expression is well-known in the 



physics literature, see, e.g., [15j. It is valid beyond the small coupling approximation [16 ]. 
We point out, however, that our formula ( fTHI) is much more general since it holds for any 
value of fi > 0. 

B. Energy Gap at Zero Temperature 

Consider now the zero temperature case T = 0. In this case, it is natural to formulate a 
functional depending only on a instead of 7 and a. In fact, for T = the optimal choice of 
■j(p) in T for given a(p) is clearly 



i(l + y/1 -4|d(p)| 2 ) for p 2 < fi 



7(P) = < • (21) 

i(l - v/l-4|a(p)| 2 ) for p 2 > /1 

Subtracting an unimportant constant, this leads to the zero temperature BCS functional 



-/ b 2 -/u| fl - a/1 - 4|d(p)| 2 ) dp+ X [ V(x)\a(x)\ 2 dx. (22) 



Defining A as in © and inserting pip , the relation between A and a at T = is 
The variational equation satisfied by a minimizer of (1221) is then 

A w=-(2^/„,^-')§i rf '' < 24 > 
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This is simply the BCS equation (jHJ) at T = 0. For a solution A, the energy gap H is defined 



as 



S = inf E(p) = inf yj (p 2 — /i) 2 + |A(p)| 2 . (25) 



It has the interpretation of an energy gap in the corresponding second-quantized BCS Hamil- 
tonian (see, e.g., [jjj] or the appendix in [lol].) 

One of the difficulties involved in evaluating H is the potential non-uniqueness of minimiz- 
ers of (|22l) . and hence non-uniqueness of solutions of the BCS gap equation (1211) . The gap 
S may depend on the choice of A in this case. For potentials V with non-positive Fourier 
transform, however, we can prove the uniqueness of A and, in addition, we are able to derive 
its precise asymptotics as A — > 0. This is the content of Lemma H] below. In particular, this 
allows us the derive an expression for H in the limit A — > 0, which is stated in Theorem [2] 

We will restrict our attention to radial potentials V with non-positive Fourier transform 
in the following. We also assume that V"(0) = (27r)~ 3 / 2 J V(x)dx < 0. It is easy to see that 

= inf spec V M < in this case, and that the (unique) eigenfunction corresponding to this 
lowest eigenvalue of is the constant function. Under these assumptions on V, we have 
the following asymptotic behavior of the energy gap 5 as A —>■ 0. 

THEOREM 2 (Energy gap). Assume that V G L 1 (R 3 )nL 3/2 (R 3 ) is radial, with V(p) < 
and V(0) < 0. Then there is a unique minimizer (up to a constant phase) of the BCS 
functional f2E) at T = 0. The corresponding energy gap, 

E = inf v / (j9 2 -^) 2 + |A(p)| 2 , 

v 

is strictly positive, and satisfies 

lim fin (£) + J[ ) = 2 - ln(8) . (26) 



a^o V \ZJ 2 y /Jib tl (X) 
Here, (A) is defined in (TT3J). 

The Theorem says that, for small A 



^ „_ e V(2VPv(A)) 

p>1 



In particular, in combination with Theorem [Tj, we obtain the following corollary of Theo- 
rem El 
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COROLLARY 1 (Universal Ratio). Under the same assumptions as in Theorem^ the 
ratio of the energy gap H and the critical temperature T c satisfies 

lim — = — « 1.7639. 

That is, the ratio of the energy gap 5 and the critical temperature T c tends to a universal 
constant as A — > 0, independently of V and /i. This property has been observed before for 



the original BCS model with rank one interaction P, llJ], and in the low density limit for 
more general interactions [15] under additional assumptions. Our analysis shows that it is 
valid in full generality at small coupling Ac 1. 

We remark that although Theorem [2] can be expected to hold under weaker assumptions 
on the potential V than the ones considered here, stronger assumptions than merely e M < 
(as in Theorem [1]) are needed for positivity of H. In particular, if A has non-zero angular 
momentum, H will, in general, vanish. 

IV. PROOF OF THEOREM Q] 

For a (not necessarily sign-definite) potential V(x) let us use the notation 

V(x) 1/2 = (sgnV(x))\V(x)\ 1/2 . 



The Birman-Schwinger principle (see Lemma 1 in [ll[) implies that the critical temperature 
T c in (flDl is determined by the fact that for this value of T the smallest eigenvalue of 

B T = XV^K^Vl 1 / 2 (27) 

equals —1. Note that although is not self-adjoint, it has a real spectrum. 

Let $ : L 1 (IR 3 ) — > L 2 (fi At ) denote the (bounded) operator that maps ip e L^R 3 ) to the 
Fourier transform of ip, restricted to the sphere Q^. Since V G L 1 (R 3 ), multiplication by 
l^l 1 / 2 is a bounded operator from L 2 (IR 3 ) to L 1 (IR 3 ), and hence 5|^| 1//2 is a bounded operator 
from L 2 (M 3 ) to L 2 (fi M ). Let 



«V( T ) = max jj— / —< \ ~~2~ j dp , 

[4irn J R3 \K T ^(p) p 2 J 



and let 

M T = K^-m^T)r$. (28) 
11 



It was shown in [ll|, Lemma 2] that V 1 I 2 Mt\V\ 1 / 2 is a Hilbert-Schmidt operator on L 2 
and its Hilbert Schmidt norm is bounded uniformly in T. (In fact, in ll| a slightly different 
definition of m M (T) was used, but it differs from ours only by a term that is uniformly 
bounded in T.) In particular, the singular part of as T — > is determined entirely by 

Since V 1 ^ 2 Mt\V\ 1 ^ 2 is uniformly bounded, we can choose A small enough such that 
sup T>0 llVVaMrlVl 1 ^!! < l/A. Then 1 + XV 1 / 2 M T \V\ 1 / 2 is invertible for any T > 0, and we 
can thus write 1 + Bt as 

1 + B T = 1 + XV 1 ' 2 {m,>(T)T$ + M T ) \V\ 1/2 (29) 

= (i + *"/»*yrn (i + 1+A ^Vi'^ 1/2rs|y|1/2 ) ' 

Then St having an eigenvalue —1 is equivalent to 



having an eigenvalue —1. The operator in (1301 is isospectral to the selfadjoint operator 

Wl 1 + \W 2 M T \V\V 2V * ' W 

acting on L 2 (r2 M ). 

At T = T c , —1 is the smallest eigenvalue of B T , hence (I30I) and (I3TI) have an eigenvalue 
— 1 for this value of T. Moreover, we can conclude that —1 is actually the smallest eigenvalue 
of (1301) and (|3~T|) in this case. For, if there were an eigenvalue less than —1, we could increase 
T and, by continuity (and the fact that m M (T) is monotone decreasing and goes to zero as 
T — > oo), find some T > T c for which there is an eigenvalue —1. Using (I2"9~j) . this would 
contradict the fact that Bt has no eigenvalue —1 for T > T c . 

Note that $V$* = y^V^ defined in ffTTl) . By assumption, = infspecV M is strictly 
negative. Since for T = T c the smallest eigenvalue of (13T]) equals —1, it follows immediately 
that 

lim Am M (T c ) = 



a^o ' inf spec $V$* v^ e M 



Together with the asymptotic behavior m M (T) ~ /i 1 / 2 ln(/^/T) as T — > 0, this imp 
the leading order behavior of ln(/i/T c ) as A — > is given by (1121) . as was proved in 



ies that 
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To obtain the next order, we employ first order perturbation theory. Since $V$* is 
compact and inf spec $V$* < by assumption, first order perturbation theory implies that 

m ^ Tc) = \(u\W$*\u) - X 2 (u\WM Tc Vr\u) + 0(A 3 ) ' (32) 
where u is the (normalized) eigenfunction corresponding to the lowest eigenvalue of $V$*. 
(In case of degeneracy, one has to the choose the u that minimizes the A 2 term in the 
denominator of (1321) among all such eigenfunctions.) 

Eq. (1321) is an implicit equation for T c . Since $VMtV$* is uniformly bounded and T c — > 
as A — > 0, we have to evaluate the limit of (u\^V M T V$*\u) as T — > 0. For this purpose, let 
(p = V$*u. Then 



(u\$VM T VT\u) 
1 



KtM 



\tp(p)\ 2 dp-m„(T) / \<p(p)\ 2 du(p) (33) 



flu 



1 [\m\ 2 - mvtp/\p\)\ 2 ] + -mvvp/\p\)\ 2 ) dp. 



KtM P 2 

Recall that K T ^(p) = \p 2 - /i\/tan\i(\p 2 - /x|/2T), which converges to \p 2 — /i| as T — > 0. 
We claim that the spherical average of |<£>(p)| 2 is Lipschitz continuous. In fact, 

(2tt) 3 I dn\^{p)\ 2 (34) 



- f dxdy Sm ^ X y ^ V(x)V(y) I e iqx u{q)duj{q) / e~ iry u{r)duj{r 



n J*? ' \p\\x-y 



and hence Lipschitz continuity in \p\ follows from the fact that V G L 1 (M 3 ), u G L 2 (Q^) and 
| sin(a)/a — sin(6)/6| < C\a — b\/(a + b) for a > 0, b > and some constant C > 0. 

Using Lipschitz continuity, it is then easy to see that one can interchange the limit and 
the radial integral over \p\, and hence obtain 

\iTn(u\WM T VT\u) = HW» , (35) 

with defined in ( fi"3l) . Moreover, this convergence is uniform in u G L 2 (fi M ). Since 
$VMtV$* is uniformly bounded in the Hilbert-Schmidt norm, this also shows that W M is 
a Hilbert-Schmidt operator, as claimed in Section IIII Al In particular, combining (132!) and 
( I35p . we have thus shown that 



mf spec [Xy/JIVf, - VW^) 
13 



It remains to calculate m M (T). Recall that K T ,n(p) = \p 2 — tanh(|p 2 — ji\/2T). We 
claim the following. 



Lemma 1. As T — > ; 



m " (T) = i^I (5^ - 7) * = ^ H + 7 - 2 + - i + ^ ) ■ (37) 



i/ere 7 7 « 0.5772 Euler's constant. 



Proof. By splitting the integral into two parts according to p 2 < /i and p 2 > /i, and changing 
variables from p 2 — fi to —t and t, respectively, we see that 

Using dominated convergence, 

For the integrals for < t < //, we use 

Jo KVJ^t ^jr+tj v * 

Moreover, using again dominated convergence, one sees that 



hm r^ZJ_y^tanh(t/2T)(ft= ^ dt = 2^ (In 2 - 1) . 

Similarly, 

lim jf* — + I - ^ tanh(t /2T) eft = 2^/1 (in 2 - 1 + V2 - In (l + V^j ) . 

In order to calculate the remaining integral t^ 1 tanh(t/2T)dt, we split the hyperbolic 
tangent into two parts, 

tanh(t/2T) = (1 - exp(-t/T)) - exp(-t/T) tanh(t/2T) . 



Using partial integration, 

p 1 — exp (-t/T) M T 1 - exp(-t) 



Jo t J 



dt 



t 



ln^ (1 -exp(-/i/T)) - / ln(t)exp(-t)dt. (39) 
T Jo 
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Since — J °° ln(i) exp(— t)dt equals Euler's constant 7 ~ 0.5772 [l7j, 4.331.1], we get 



Finally, we have to evaluate 

tanh(t/2T) /"^ tanh(t/2) 



tanh(t/2T) , /m , , tanh(t/2) , , , 
J Y -exp(-t/T)dt = j p—*- exp(-t) dt . 



In the limit T — > 0, this becomes 



13, 3.411.28] 

p tanh(t/2) . , , , 7T 

/ ^— ^exp(-t)dt = In-. 

Jot 2 

Collecting all the terms, we arrive at ( 1371) . □ 



Theorem [T] follows immediately from fl36|) and fl37j) . recalling the definition of ^(A) in 
(ED and ([15]). 

V. PROOF OF PROPOSITION [Tj 

For w(p) = (47r / u)~ 1 / 2 the (normalized) constant function in L 2 (f2 At ), we have 

|/" (V(p-q) ^^\MP)MQ) 

x 2 r / 1 1 



(u|B» - a (A) = A / (>(p - q) - V(0) 

V ^ Jn u xn u v 



4vr J R 3 VlP 2 -/^l P 2 



47T/i 47TyU 

(P)\ 2 ~ \M^P/\P\)\ 2 ) dp 



+ ~ I -2 (\MP)\ 2 - \MP)\ 2 ) dp , (40) 



^ [ 1 

4vr 7 R 3 p 2 

where V M (p) = V"(p - q)du IM (q) / (Anfi) and Vo(p) = ^(p). 
Consider first the term linear in A. It is given by 

,2/ 



Air 



For < /3 < 2, ( v //I|x|)~ /3 (sin 2 (y / /i|x|)/(/i|a;| 2 ) — 1) is a function that is bounded indepen- 
dently of fi, and goes to zero as \i — > for every x. Hence, if J \ V(x)\\xfdx < 00, it follows 
from dominated convergence that 

[ V{x) ( Sin2{ ^ Xl) -l)dx = 0. 
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The second term in (HOl) can be rewritten as 



A 2 / v ^ Mvmx\) v ^MVfj\y\) r ( p 2 



v^M y/i*\v\ Jo vb 2 -/^ 

x / sin(p|x|) _ sin(^|x|)\ / sin(pM) + sin(^\y\)\ 

V p\ x \ vWl J \ p\y\ vwl J 

The last term in the integrand is bounded by 2, while the second term is bounded by 
C\p — ^JJi\ min{|x|, l/(p + y/Ji)}, since | sin(a)/a — sin(6)/6| < C\a — b\ min{l, 1/(6 + a)}. 
Since < (3 < 2, we can estimate min{|x|, l/(p + y/Jt)} < \x\@/ 2 (p + y^u) _1 - Hence we 
conclude that (1411) is bounded by 

P 2 ,\ \P~\M 



Since 1 + [3/2 > [3 this yields a bound of the desired form. 
Finally, the last term in ( 1401) is given by 

V 47r / Jr« F-y| V v^fI v^M / 

The Hardy-Littlewood-Sobolev inequality [18| and Holder's inequality imply that 

/ V ,^y kl^ /2 lyl^dy < C|[y| • f% b <C\\Vh /2 J \V(x)\\x\Pdx 

for some positive constant C. Moreover, since for < (3 < 2 the expression 
|y|) _/3 / 2 (l — sin(y^|x|) sin( v //I|2/|)/(/ x l a; ll2/|)) is uniformly bounded and goes to zero 
pointwise as fi — > 0, it follows again from dominated convergence that 

lim J_ r V(x)V(y) ( snuVgN) sin(^) \ _ Q 
i^* y//fJw> \x-y\ V V^kl v^M / 

This proves Proposition [TJ 



VI. PROOF OF THEOREM [2] 

We start by showing that the minimizer of the BCS functional (1221) at T = is unique 
under the assumption that is non-positive and V(0) < 0. Note that if, in addition, V is 
radial, this necessarily implies that also the minimizer has to be radial. 
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Lemma 2. Let V G L 1 n L 3 / 2 and assume that V < and V(0) < 0. Then the BCS 
functional JF at T = 0, defined in (CHjj, has a unique minimizer a (modulo a constant 
phase), whose Fourier transform a is strictly positive. 

In particular, this implies that A is strictly positive. From the gap equation (I24p one 
easily concludes that A is also continuous, and hence S > 0. 



Proof. The existence of a minimizer was shown in [10|. Moreover, since e M = inf spec V M < 
under our assumptions on V, the critical temperature T c is strictly positive, and hence a 
minimizer of T§ is necessarily not identically zero. 
Note that since V < 0, 



Hence, if a(p) is a minimizer of JF , so is |a(j>)|. 

Assume now that there are two different minimizers / ^ g, both with nonnegative Fourier 
transform. Since t i— ► 1 — \/l — 4t is strictly convex for < t < 1/2 we see that ip = 



This is a contradiction to /, g being distinct minimizers, and hence / = g. In particular, the 
absolute value of a minimizer is unique. 

Let a be the unique non-negative minimizer. The corresponding A in ( 1231) is also non- 
negative, and satisfies the gap equation ff24l . From this equation it follows easily that A 



is strictly negative in a non-empty open ball around the origin. From (124"]) it follows that 
A can only vanish at a point po if it vanishes in this open ball centered at po, and hence is 
identically zero. 

In particular, from (123]) we conclude that also a is strictly positive. Since we already 
know that all minimizers must have the same absolute value, this implies that any (not 
necessarily positive) minimizer of Tq is non-vanishing. But the inequality (j4"2"l) is strict for 
non-vanishing functions, unless &{p) = e lK '\a(p)\ for some constant n£l. This proves the 
uniqueness of the minimizer. □ 

In the following, we shall choose the arbitrary constant phase factor in the (otherwise) 
unique minimizer of JF as 1, i.e., we take a to be positive. 





(42) 



■j=f + i-j=g, satisfies 



W) < IMf) + \?o{9) ■ 



is strictly positive. In fact, since V(p) is continuous and strictly negative at the origin, it 
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Note that the variational equation for the minimizer of T§ can be written as 

(E(-iV) + \V(x))a(x) = 0, (43) 

with E(p) = \p 2 - fi\/y/l -4\a{p)\ 2 = ^\p 2 - fi\ 2 + \A{p)\ 2 , see p3l and (J23j| . That is, a is 
an eigenfunction of the pseudo differential operator E{— iV) + AK(x), with zero eigenvalue. 
Under our assumptions on V, we can even conclude that it is the ground state. 

Lemma 3. Let V < and V(0) < 0. Lei a be the minimizer of the BCS functional |H)) ; 
with corresponding A defined in Then a is the ground state of the operator 

E(-iV) + XV(x) , (44) 

where E(p) = a/ (p 2 — fi) 2 + \A(p)\ 2 . 

In particular, this implies that E{— iV) + XV(x) > 0. This property will be an essential 
ingredient in the proof of Theorem [2j 

Proof. Since V(p) < 0, the ground state of (1441) can be chosen to have non-negative Fourier 
transform. It is therefore not orthogonal to a, since a is strictly positive by LemmaEJ Hence 
a must be a ground state. □ 

Similarly to the proof of Theorem [1] in Section IIV1 we can employ the Birman-Schwinger 
principle to conclude from Lemma [3] that <p — V 1 ^ 2 a satisfies the eigenvalue equation 

XV 1 ' 2 - 1 |V| 1 / 2 0=-0, (45) 



Moreover, there are no eigenvalues smaller than —1 of the operator on the left side of (j45j) . 
Let 

fhJA) = max I - — / [ ; - — ) dp , I . (46) 



Similarly to ( 1281) . we split the operator in ( 1451) as 
E(—i\) 



By proceeding in the same way as in the proof of 
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Lemma 2], one shows that V 1 / 2 M/^\V\ 1 / 2 



is bounded in the Hilbert-Schmidt norm, independently of A. Moreover, as in the proof of 
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Theorem □ (cf. Eqs. fl22H2H)), the fact that the lowest eigenvalue of \V 1 l 2 E(-iVy 1 \V\ 1 / 2 
is —1 is, for small enough A, equivalent to the fact that the selfadjoint operator on L 2 (f2 M ) 

y' /! m^ y ' <i8, < 47 > 

has —1 as its smallest eigenvalue. 

Recall that $V8* equals ^[JiV^ defined in (TTTT) . Our assumptions on V imply that 
the lowest eigenvalue e M of V M is strictly negative, and non-degenerate. This implies that 
lim^o Am M (A) = — 1/( V / /Ie /J ) and hence, in particular, m M (A) ~ A -1 as A — ► 0. The 
unique eigenfunction corresponding to the lowest eigenvalue < of is, in fact, a 
positive function, and because of radial symmetry of V it is actually the constant function 
u{p) = (4 7 r / u)- 1 / 2 . 

We now give a precise characterization of A(p) for small A. 

Lemma 4. Let V E L 1 D L 3 ^ 2 be radial, with V < and V(Q) < 0, and let A be given in 
&3\) . with a the unique minimizer of the BCS functional (TJS|). Then 



A(p) = -/(A) / l/(p-g)da;(g)+A^(p) (48) 



for some positive function f{\), with ||^a|U°°(]r 3 ) bounded independently of X. 

Proof. Because of (]4"5l) . t?| ^| 1 ^ 2 <^> is the eigenfunction of (1471) corresponding to the lowest 
eigenvalue —1. Note that because of radial symmetry, the constant function u(j>) = (Att/j,)^ 1 ^ 2 
is an eigenfunction of (j4"Tj) . For small enough A it has to be an eigenfunction corresponding 
to the lowest eigenvalue (since it is the unique ground state of the compact operator $V$*). 
We conclude that 

41 = f{X \ + XV^M A \V^ V '' 2rU = /(A) < V " /2r " + A& > (49) 
for some normalization constant /(A). Note that ||£a||2 uniformly bounded for small A, since 
both V 1/2 M A \V\ 1/2 and V 1/2 ^* are bounded operators. 

Recall from ([23]) and and the definition = V 1/2 a that 

A(p) = 2E(p)a(p) = -2XVa(p) = -2\\V\V*<j>(p) . 
In combination with (j49|) this implies that 

A(p) = -2A/(A) (v5^(p) + \vx(p) 
19 



with rjx = IV] 1 / 2 ^. Note that V$*u(p) equals (47r/i) -1 / 2 J Q V(p — q)dio(q). Moreover, since 
II^aIIoo < (27t)~ 3 ^ 2 ||^a||i < (27I")~ 3/ ' 2 ||V a ||i||^a||2 by Schwarz's inequality, ||^a||oo is bounded 
uniformly in A. This implies the statement of the Lemma. □ 

Note that J Q V(p — q) du(q) is a Lipschitz continuous function. In fact, 

/ (v{p-q)-V{r-q)]du{q) = \l^ ! V(x) Sm( /^ |x|) ( e ~^ - e ~ irx ) dx 
Jn„ v > V 7T J R3 \x\ 

< \ —\p — r \ [ \V(x)\\sxn(^/u\x\)\dx . 
V vr J R3 

Using this property, together with Lemma HI we can now estimate rh^(A) in fT46]) . We first 
consider A to be a constant, and start with the following observation. 

Lemma 5. Let ■d > 0. As ■d — > ; 

m (7?) = - — / ( . ? I <fo = — (ln^-2 + ln8 + o(l)) . (50) 

Proof. Proceeding as in the proof of Lemma [U we have 



"m(*) = ^/ 1 V " m m V - -^-^)^ (51) 



The first integral becomes (In 4 — 2 — ln(l + y/2))/ \fii in the limit $ — > 0. The last integral 
becomes /i - 1 / 2 ln(l + V2). Finally, 



= v^ln = ^ln — + o(l) . 

This proves the statement. □ 

Since m M (A) — > oo as A — > 0, if follows from Lemma H] that Km^o f(ty = 0. In the 
following, we shall slightly abuse the notation and denote by A(y / /T) the value of A(p) on 
the Fermi sphere f2 M . Since A is a radial function, this is well defined. Using Lemma [5], we 
will now argue that 

fi " (A) = ^( ln A(^)- 2 + " l8 + 0(1) ) (W) 
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as A — > 0. In fact, by inspection of the proof of Lemma [5] we see that all we have to show is 
that 

-dt+- / dt = In -— -— — - + o(l) , 



which follows easily from Lemma HI Lipschitz continuity of j n V(p — q) du(q) and the fact 
that liniA-^o /(A) = 0. 

From (1471) we conclude that 

m ^ (A) = \(u\W$*\u) - \ 2 (u\WM A V$*\u) + 0(A 3 ) ' (53) 

where u(p) = (47r/i)~ 1//2 is the normalized constant function on the sphere f2 M . Moreover, 
with ip = V$*u, 

(u\wm a vt\u}= I -^—i^dp-m^A) [ m^p/\p\)\ 2 Mp) 

f 1 r, s,o , _ on 1 



[W) [mp){2 ~ '^ p/|p|)|2] + pMv^p/\p\)\ 2 ) dp. 

Using Lemma H] and the fact that lim^o / (A) = 0, we conclude that 

\im(u\WM A V$*\u) = (u\W„\u) , (54) 

A — *0 

with W M defined in (ITS]) . (Compare with Eqs. (1531) and (1551) .) In combination with (I5"2"j) and 
(I55|) and the definition of 23^ in (fT4l) . this proves that 

lim fin ( -r^—) + 7^ , r l = 2 - ln(8) . 

The same holds true with (u\Bfj,\u) replaced by 6 M (A) = inf spec £> M , since under our assump- 
tions on V the two quantities differ only by terms of order A 3 , as explained in Section fill Al 
Now, by the definition of the energy gap 5 in (T25]) . 5 < A(^//J). Moreover, 

S > min I A(p) \ , 

b 2 -H<= 

from which it easily follows that H > A( v //I)(l — o(l)), using Lemma HI This proves 
Theorem [2j 
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VII. CONCLUSION 



We have presented a rigorous analysis of the critical temperature T c and the energy gap 
5 for the BCS gap equation. Our results are valid for a general class of local interaction 
potentials, and for any value of the chemical potential. The correct expressions for T c and 
S, valid to second order in the Born approximation, do not seem to have appeared in the 
literature before. 
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